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Results are presented for the nonequilibrium dynamics of a quantum XXZ-spm chain whose 
spins are initially arranged in a domain wall profile via the application of a magnetic field in the 
z-direction which is spatially varying along the chain. The system is driven out of equilibrium in 
two ways: a), by rapidly turning off the magnetic field, b). by rapidly quenching the interactions 
at the same time as the magnetic field is turned off. The time-evolution of the domain wall profile 
as well as various two-point spin correlation functions are studied by the exact solution of the 
fermionic problem for the XX chain and via a bosonization approach and a mean-field approach 
for the XXZ chain. At long times the magnetization is found to equilibrate (reach the ground 
state value), while the two-point correlation functions in general do not. In particular, for quenches 
within the gap less XX phase, the spin correlation function transverse to the z direction acquires 
a spatially inhomogeneous structure at long times whose details depend on the initial domain wall 
profile. The spatial inhomogeneity is also recovered for the case of classical spins initially arranged 
in a domain wall profile and shows that the inhomogeneities arise due to the dephasing of transverse 
spin components as the domain wall broadens. A generalized Gibbs ensemble approach is found to 
be inadequate in capturing this spatially inhomogeneous state. 

PACS numbers: 71.10.Pm, 02.30.1k, 37.10.Jk, 05.70.Ln 



I. INTRODUCTION 

Recent years have seen remarkable progress in manipu- 
lating cold-atom systems^ providing us with almost ideal 
realizations of strongly correlated many-particle systems. 
Cold atom systems also have the unique feature that the 
interactions between particles and the external poten- 
tials that they are subjected to are highly tunable and 
can be changed rapidly in time, thus driving these sys- 
tems into highly nonequilibrium states^. This property 
has in turn motivated a lot of theoretical activity that 
revolves around studying nonequilibrium time-evolution 
of strongly correlated systems arising due to a sudden 
change in some parameter of the Hamiltonian referred to 
as a "quantum quench" . The general consensus is that 
if the strongly correlated system is integrable, its time- 
evolution from some arbitrary initial state is highly con- 
strained by the initial conditions so that these systems 
do not reach the ground state but instead reach interest- 
ing nonequilibrium time-dependent or time-independent 
steady states. 

In this paper we study quenched dynamics in one such 
integrable model, namely the XXZ spin chain. This 
model shows rich behavior even in its ground stated. For 
example, it exhibits a quantum critical point at J z = ± J 
where J is the exchange interactions for the x, y com- 
ponents of the spins, and J z is the exchange interaction 
for the z component of the spins. For \J/J Z \ > 1, the 
ground state is an XX phase characterized by gapless 
excitations, while for \J/J Z \ < 1, the spins are in an 
Ising phase where the excitation spectrum has a gap. 
Note that the special point in the gapless phase corre- 
sponding to J z = will be referred to as the XX-model, 
while the term gapless XX phase will refer to the regime 



\J/J Z \ > 1. 

Nonequilibrium dynamics of the XXZ model arising 
due to quenches has been a very active area of research 
(see for example^ ). While most previous work has stud- 
ied quenched dynamics in spatially homogeneous sys- 
tems, in this work we consider the nonequilibrium time- 
evolution arising from quenching from an initial state 
which is spatially inhomogeneous. Such an initial state is 
created by the application of an external magnetic field in 
the z direction that varies in magnitude along the chain, 
changing its sign at some point so that the spins are 
aligned in a domain wall pattern. (SeeS for an experimen- 
tal realization of such a setup) . We study nonequilibrium 
dynamics that arise in the following two ways: a), by a 
sudden quench of the magnetic field to zero; b). by a 
sudden quench of the magnetic field to zero along with 
a quench in the magnitude of the exchange interactions. 
Note that the above are situations where the time evo- 
lution is due to a final Hamiltonian that has a spatially 
homogeneous ground state. We would like to understand 
how the system evolves in time after the quench, and 
whether the initial inhomogeneity affects the properties 
of the system at long times. 

Some literature already exists which studies time evo- 
lution of a spin chain initially arranged in a domain wall 
profile. Antal et ali study the time evolution of a sharp 
domain wall in the exactly solvable XX model and find 
ballistic broadening of the domain wall width. Simi- 
lar studies for the time-evolution of a domain wall has 
been carried out for the XX model^— and the critical 
transverse spin Ising mode&ii. The numerical method 
time-dependent density matrix renormalization group (t- 
DMRG) has been used to study the time-evolution of a 
domain wall both in the gapless phase as well as the 
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gapped Ising phase of an XXZ chai n 12 i 13 . Their study 
reveals qualitatively different behavior in the two phases, 
with ballistic broadening of the domain wall in the gap- 
less phase, and more complicated non-ballistic behav- 
ior in the Ising phase. In addition conformal field the- 
ory methods^ have been used to study both the time- 
evolution of the domain wall profile as well as two-point 
correlation functions for the Ising chain. Studies also ex- 
ist on the time-evolution of other quantities such as the 
entanglement entropy after a quench from a spatially in- 
homogeneous initial stat o 15 i 16 . Thus what is lacking in 
the literature is the study of the time evolution of two- 
point correlation functions for the XXZ chain from an 
initial inhomogeneous state. The advantage of studying 
these quantities is that they are often more sensitive to 
nonequilibrium initial conditions than averages of local 
quantities such as the magnetization. 

In this paper we extend previous results by studying 
the time evolution of both the local magnetization as well 
as two point correlation functions related to the longitu- 
dinal (along z direction) and the transverse spin corre- 
lation function for the XXZ chain. In addition we use 
a variety of theoretical methods such as bosonization, 
exact solution of the fermionic problem as well as a clas- 
sical treatment of spins which allows us to compare the 
relative merits of these approximations. Consistent with 
previous results, we find that the magnetization always 
equilibrates (i.e, reaches the ground state value). In ad- 
dition the details of how the domain wall profile evolves 
in time depend on whether the system after the quench is 
in the gapless phase or the gapped Ising phas o 12 ' 13 . For 
example in the gapless phase the domain wall is found to 
spread out ballistically, however in the gapped phase the 
dynamics is much more complicated showing oscillations 
and revivals at short time scales. 

In studying the time-evolution of two-point correlation 
functions we obtain the interesting new result that the 
transverse spin-correlation function shows a lack of equi- 
libration reaching a nonequilibrium steady state that is 
also spatially inhomogeneous. The length scale of the 
spatial inhomogeneity is found to depend on the details 
of the initial domain wall profile and is found to arise due 
to the dephasing of transverse spin components as the do- 
main wall broadens. This result is qualitatively different 
from those obtained for the Ising chain initially arranged 
in a domain wall profile where no residual inhomogeneity 
was observed^. 

The paper is organized as follows. In section [TT] we 
study the XX model which is initially subjected to a 
spatially varying magnetic field so that the spins are ar- 
ranged in a domain wall profile. We study the nonequi- 
librium dynamics that arises when this magnetic-field is 
suddenly switched off. These results are obtained for 
two cases. One for an initial external magnetic field that 
varies linearly in position and is equivalent to an effec- 
tive electric field on the fermions. The second case is 
that of an initial magnetic field which has a sharp step- 
function profile. In section IIIII the effect of quenching 



only a spatially varying magnetic field in the XX model 
is studied via a bosonization approach. The results of 
section [TT] and UTTJ are found to be in qualitative agree- 
ment and in particular recover ballistic domain wall mo- 
tion and spatial inhomogeneities in the transverse spin 
correlation function. A physical explanation for this in- 
homogeneity is provided in section lTVl where the problem 
of quenching a spatially varying magnetic field is studied 
classically and a spin- wave pattern is obtained as a long 
time solution of a Landau-Lifshitz equation. 

In section |Vl the effect of simultaneously quenching 
the magnetic field and interactions is studied using a 
bosonization approach. In section IV Al results are pre- 
sented for the case when the quench is entirely within 
the gapless XX phase. The spatial inhomogeneities in 
the correlation functions are recovered along with the 
nonequilibrium exponents obtained earlier in a purely ho- 
mogeneous interaction quenc h 17 i 18 . In section [VB] the ef- 
fect of quenching into the gapped Ising phase is studied 
via a semiclassical approach. Here we find the domain 
wall dynamics to be qualitatively different, however at 
long times we find that all inhomogeneities eventually 
decay away. The latter result may very well be an arti- 
fact of the semiclassical approach which neglects creation 
of solitons. The results of the semiclassical approach is 
complemented by a mean-field treatment in section IVII 
where we address the question of whether Ising order 
can develop for a system initially in a domain wall state. 
Finally we conclude in section IVIII 



II. QUENCH OF A SPATIALLY VARYING 
MAGNETIC FIELD IN THE XX MODEL: EXACT 
SOLUTION OF THE FERMIONIC PROBLEM 

In this section we will study the XX model with a 
magnetic field in the z direction that varies linearly in 
position and changes sign at the center of the spin chain, 
so that the ground state is a domain wall configuration. 
In subsection III Al we will study the properties of this 
domain wall state, and in subsequent subsections study 
the nonequilibrium time-evolution arising when the the 
spatially varying magnetic field is suddenly switched off. 



A. Creation of domain-wall state 

Our starting point is the XX Hamiltonian in a mag- 
netic field, 
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where Sj — a v are Pauli matrices, hj is an ex- 

ternal magnetic field aligned along the z direction, and 
J is the exchange energy. For a uniform magnetic field 
one commonly employs the Jordan- Wigner transforma- 
tion-!^— to map the spin system to a system of spin-less 
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non-interacting fermions which can be easily diagonalized 
(for details see Appendix [3}. 

To create an inhomogeneous state we will consider a 
linearly varying field hj = Fja, where a is a lattice spac- 
ing. Note that the ground state properties of spin-chains 
in various spatially varying magnetic fields has been stud- 
ied ir*22. After the Jordan- Wigner transformation Eq. Q] 
for hj — Fja maps to the Wannier-Stark proble m 15 i 24 i 25 
of electrons in a lattice subjected to a constant force (or 
electric field) J 7 , 



Cj+l 



c j+i c j 



ac'jCj (2) 



This Hamiltonian can be diagonalized as follows, 

H X x — ^ ^rnVm'Tlmi 



(3) 
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mJ-a, (V m = integer) 



(4) 



FIG. 1: (Color online) Plot of the initial magnetization for 
various values of the domain wall width a = S- . 
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Above Cj — -jjj ^2 k Cke lk \ N is the total number of sites 
and J n (x) is a Bessel function of the first kind. 

We construct the ground state |$) by including all neg- 
ative energy states, 



i*) = n ^ io> 



(7) 



m<0 



This is a domain-wall state with a characteristic width 
- = a = to the left (right) of which all spins are up 
(down). This can be seen by evaluating the magnetiza- 
tion at the position ja 



m z {j) = h\c] Cj \A - ~ 



(8) 



which on rewriting Cj , cj in terms of the ry m operators 



and evaluating the expectation value is found to be 

™ 2 w = -o + > :^+m(«) (9) 



m>0 



Upon noting that the support of J n (x) is restricted to2S 
\n\ < x, and using the identity 2 ^ J2 n ^n( x ) = U ^ is eas y 
to see that this indeed describes a domain-wall state of 
width a. Eq.[9]is also plotted in Fig.Q] It is worth noting 
that the case of a = corresponds to an infinite "electric 
field" T that forces a sharp wall of zero width. This case 
was investigated by Antal et a£. We will extend some 



of their results to more general domain walls, and also 
study spin-spin correlation functions. 

The basic spin correlation functions that we will com- 
pute are, 

C» v (j,j + n) = (S%S» +n ), {v = x,z) (10) 

where we will refer to C zz (C xx ) as the longitudinal 
(transverse) spin correlation function. These may be ex- 
pressed in terms of Majorana operators^ Aj = cj + Cj, 
and Bj — cj — Cj as follows, 

C zz (j,j +n) = i (B 3 A 3 B J+n A 3+n ) , (11) 

and 

C xx {jJ+n) = - (BjAj +1 B j+1 ■ ■ ■ Aj +n -iBj +n -iAj +n ) , 

The above correlations can be evaluated by rewriting Ci 
in terms of r] m and applying Wick's theorem. The basic 



contractions are 



2,s 



(BjAj +n ) — - (AjB j+n ) — N [J 3 + m (a>)J 3+n+m (a) 



m>0 



Jj— m\P^)Jj-\-n— m(^v] i 

(AjAj +n ) = — (BjBj +n ) = S n=0 . 



(13) 
(14) 



Using some Bessel function identities, the mixed contrac- 
tion can be simplified to 

a 

(BjAj+n) = [ J j+n{a)J 3+ i{a) - J j+n+1 {a)J 3 (a) 



-Jj +n (a)Jj-i(a) + Jj +n -i(a)Jj(a)} . 



(15) 



Since we encounter no contractions of the form (AiAi), 
the function C xx (j, j + n) may be written^ as a Toeplitz 
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FIG. 2: (Color online) Plot of C zz {j,j + n) for n = 1,2,3,4 FIG. 3: (Color online) Plot of C xx (j,j + n) for n = 1,2,3,4 



and a = S- = 50. 



and a = S- = 50. 

J- a 



determinant which is a computationally cheap task. In 
the next sub-section when we evaluate the correlation 
functions after a quench of the magnetic field, we will 
find (AjAj +n ) =/= in the nonequilibrium state. In this 
case, the evaluation of C xx will require computing the 
square root of a determinant^!. 

The correlation functions C zz and C xx for a = 50 
are given in Fig. [5] and Fig. [3] respectively. The main 
point to notice here is that close to the center of the do- 
main wall where the external magnetic field is small, the 
correlations tend to mimic those in the ground state of 
the homogeneous (hj = 0) XX model. For example, the 
nearest neighbor longitudinal spin correlation function in 
the absence of a magnetic field is22 C zz (n = 1) = — 
This is precisely the value that the C zz (j,j + 1) correla- 
tion function takes in Fig. [3] at the center of the domain 
wall. 



Domain wall dynamics after the magnetic field 
quench 



The time-evolution of Cj , ct in terms of the r\ m operators 
is given in Eqns. 1X61 IA7I from which we obtain 



'(j,t) = -\+ J2 \W + m,t,a)\ 2 , 



m>0 



where 15 

L(j + m,t,a) 



n dk 



(18) 



(19) 



Evaluating the integral, we find 

m'(j,t) = 4+E J l+m (VW+a 2 ) ■ (20) 

m>0 

T his corresp onds to a wall whose width W = 
(Jt) 2 + a 2 increases linearly and hence ballistically in 
time with a velocity of J a for times t > a/ J. Setting 
a = 0, we recover the result of Antal et a£. Note that 
the time-evolution of the entanglement for general a has 
been studied ir>i£. 



We will now explore the dynamics that arises when 
the spatially varying magnetic field is suddenly switched 
off at t = i.e., hj{t) = 6(-t)Fja. Thus at t < the 
many-body wavefunction of the system is Eq. [71 while 
for t > 0, as shown in Appendix [XI the wave-function 
evolves according to the Hamiltonian, 



(16) 



where = — Jcos(fca). 

The time-dependent magnetization is given by 



m z (j,t) 



<*|ct(t) Ci (t)|*) 



(17) 



C. Correlation functions after the magnetic field 
quench 

We now turn to the evaluation of the longitudinal and 
transverse correlation functions given in Eqns. f)l 1 1) and 
(TT21 at a time t after the quench. The basic contractions 
that we need are (Bj(t)Aj +n (t)) and (Aj(t)Aj +n (t)) for 
which we find the following expressions 

(Bj(t)A j+n (t)) = 

^ [e mS (J J+n (r)J J+1 (r) - J j+n+1 (r)Jj(r)) 
-e~ mS {J 3+n {r)J 3 ^{r) - J ]+n ^{r)J 3 {r))} , (21) 
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(A 3 (t)A 3+n (t)) 



e mS (J j+n (r)J j+1 (r) - J i+n+1 (r)J,(r)) 



2n 



+e~ inS (J 3+n (r)J 3 ^(r) J 3+n ^{r)J {r))] , (22) 



with r = y/ (Jt) 2 + a 2 and tan 6 = Jt/a. As be- 
fore, we have (B 3 (t)B 3+n (t)) = — (Aj(t)Aj +n (t)) and 
(B j (t)A j+n (t))=-(A j (t)B j+n (t)). 

The transverse correlation function C xx (j,j + n,t) is 
given by 



1 



\C xx (j,j+n,t)\ = -VMC, 
where C is the anti-symmetric matrix 

with the n x n sub-matrices defined by 



(B 



i+j 



-iB. 



k+j-i/ 



Sik — Ski — 

Qik = -Qki = (A l+j A k+3 ) (k > i) 

Gik = \Bi+j-lAk+j) . 



(k > i) 



(23) 



(24) 



(25) 
(26) 
(27) 



Note that both S and Q are anti-symmetric, and thus 
defined by the elements with k > i, while no such restric- 
tion is placed on G. 

The numerical evaluation of the above expression for 
a — 25 is shown in Fig.[4]for three different times: t = 0, 
an intermediate time and for long times where one finds 
the appearance of a spatially oscillating pattern at the 
scale of the Fermi- wave vector kp — ir/2. 

In order to get some insight into the long time behav- 
ior of the system, we consider the limit of r — > oo while 
imposing the condition j, j+n <ron the two position in- 
dices of the correlation functions. This limit corresponds 
to a very broad domain wall which may arise either due 
to a very weak electric field a > 1 at arbitrary times 
t (and hence arbitrary S) or could arise on waiting long 
enough after a quench so that an initial narrow domain 
wall has had enough time to broaden. The condition on 
the position indices imply that we are looking at corre- 
lations in the vicinity of the center of the domain wall 
where the magnetization has locally equilibrated. Using 
the asymptotic expansion for the Bessel functions with 
large arguments^, we find 

/t, r x 4 , %, 2cos(n<5) . /n7r\ 
{B 3 (t)A ]+n (t)) ±-L sin — , (28) 

TTTL V Z / 

2i sin(7i(5) 



{A 3 {t)A 3+n {t)) 



■ sm 



/n7T\ 



(29) 



Substituting the asymptotic expansions in Eq. 
Eq.[29]we find, 



C**(j,j + n,t)^CI*(n) = -'- 



V 2 ) 



(30) 



where C^(n) is the equilibrium result (i.e. the result in 
the ground state of Hf) for the ab correlation function 22 . 
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FIG. 4: (Color online) Time evolution of C xx (j = 0, n) cor- 
relation function for a = 25. Solid line is the equilibrium 
correlation function Ctl ' . 



Employing the asymptotic expansions in the evalua- 
tion of C xx (j, j + n), we have verified through n ~ 1000 
that 



C xx (j,j+n,t) =C xx (n)cos(nS) 



eq 
t—s-oo 



¥ C xx (n) cos 



/irn\ 



(31) 
(32) 



where C xx (n) is the equilibrium value of the xx- 
correlation function 19 which has the asymptotic form 29 , 



CTM 



with logK = i/ °° f (e" 



(33) 



cosh 2 t ) 



The appearance of these oscillations at a wave-length 
of A = 4 where 

C xx (. h j + n, t) ^> C xx (n) cos(^) (34) 

is indeed intriguing. To understand better what sets this 
length scale, we study more general initial domain walls. 
The one considered so far has the maximum possible pos- 
itive (negative) polarization of S z — ±1/2 on the left 
(right) ends of the domain wall, with a width of the do- 
main wall controlled by a — . Antal et at have shown 
how to construct domain walls of zero width (a = 0) but 
arbitrary polarization ±mn at the two ends. Starting 
with this initial state i.e. 



S z At = 0) = -m sgn(j) 



(35) 



we study the time-evolution under the XX model and 
extend the results of Antal et al to the study of the trans- 
verse spin correlation function. 

Details of the computation are given in Appendix |5] 
In the long time limit of Jt — ¥ oo, analytic expressions for 
contractions in Eq . |2"51 |2T)1 generalize to Eq. lB9l These can 
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then be used to construct the determinant of C (Eq. 
required for the computation of C xx . The results for 
C xx in the long time limit for several different domain 
wall heights too are shown in Fig. [5] Indeed we see the 
appearance of a spatial oscillation where the wavelength 
depends on the height too of the domain wall as follows, 



A = — 

too 



(36) 



Thus the main results of this section on the effect 
of quenching a spatially varying magnetic field in a 
XX model are as follows: a). An initial domain wall 
spreads out ballistically after the magnetic field quench 
so that the magnetization equilibrates to its homoge- 
neous ground state value of m z — 0. For a sharp do- 
main wall, (a = y-^ -C 1), the magnetization at a dis- 
tance na from the center of the domain wall begins to 
equilibrate after a time Jt ~ na. This is an example 
of the horizon effec t 30 i 31 , namely the minimal time re- 
quired for ballistically propagating excitations from the 
center of the domain wall to reach the observed point; 
b). Another consequence of the horizon effect is that 
any two point correlation function C ab (j = 0, n) is found 
to change significantly only after times 3 -! Jt > \n\ for 
a <C 1; c). The C zz (j,j + n) correlation function even- 
tually equilibrates (i.e. reaches the ground state value) 
at long times where the approach to equilibrium is a 
power-law of C(-^); d). The C xx (j,j + n) correlation 
function reaches a nonequilibrium steady state at times 
min(\J a 2 + (Jt) 2 , fJa) > 1 to a value given by Eq. [34] 
which is basically the ground state correlation function 
with oscillations at a wave-length A = ^- superimposed 
on it, where to is the height of the domain wall, i.e., the 
magnitude of the maximum polarization at its two ends. 

Note that the effect of quenching from an initial do- 
main wall state was also studied in^ for the critical trans- 
verse spin Ising model (J2i SfS x +1 + J2i S() using the 
methods of boundary conformal field theory. Our results 
for the XX model differ from those ofi^ in two impor- 
tant ways. One is that at long times we recover a spatial 
spin-wave pattern that was not captured ini^. Secondly 
at long times after the quench, our two-point correlation 
functions continue to be critical (i.e., have a power-law 
decay) whereas the results o£i^ point to thermal behav- 
ior with exponential decay in the two-point correlation 
function. The differences between these results are due 
to the very different universality classes and conserva- 
tion laws of the XX model and the transverse spin Ising 
model. For example, while the total magnetization in the 
z-direction is conserved in the XX model, it is not con- 
served in the Ising model. As a result the magnetization 
has a nontrivial time evolution in the Ising model even 
for points outside the domain wall horizon, which even- 
tually gives rise to a thermal behavior. In addition, as 
we shall show in section IIV1 the appearance of the spin- 
wave pattern in the transverse spin correlation functions 
in the XX model is due to precession of spins in an ef- 
fective magnetic field created by the initial domain wall. 



This precessional physics is absent for Ising spins. 

We would also like to point out that the spatial os- 
cillations found by us have also been observed in a dif- 
ferent physical setting by Rigol et a£2 who have studied 
inhomogeneous quantum quenches for one-dimensional 
hard-core bosons. Via a Jordan- Wigner transformation 
hard-core bosons in one dimension can be mapped onto a 
system of spinless fermions. The physical situation con- 
sidered by Rigol et al was one where an initial strong con- 
fining potential produces a local region of high density. 
(In the spin language this would correspond to creating 
a pile up of magnetization in the center of a spin-chain). 
The time evolution of the system when this confining 
potential was suddenly switched off was studied. The 
initial inhomogeneity was found to move out ballistically 
in two opposite directions. Moreover, the analog of the 
transverse spin correlation function which in the boson 
language corresponds to the off-diagonal matrix element 
was found to show spatial oscillations in the two out- 
going lobes. (This is somewhat different from our set-up 
where the spatial oscillations appear at the center of the 
chain and extend over longer portions of the chain as the 
domain wall broadens). For an initial tight confinement, 
Rigol et al found that the oscillations in the lobes ap- 
peared at precisely the same wave-length of A = 4 that 
we find here for a domain wall of height too — 1/2. For 
the case of bosons, these oscillations had the physical in- 
terpretation of the appearance of quasi-condensates at 
the wave-vector Icf = ±7r/2. A similar observation that 
an initial spatial inhomogeneity for a system of hard core 
bosons can lead to correlations in momentum space was 
also made by Gangardt et a^ 3 .. 

In subsequent sections we will study the effect of 
quenching a spatially varying magnetic field with and 
without an interaction quench in a XXZ chain using a 
bosonization approach. For quenches within the gapless 
XX phase, results of this section, namely critical behav- 
ior in the two-point correlation function even after the 
quench, and the presence of spatial inhomogeneities will 
be recovered. In addition we will find that the wavelength 
of the inhomogeneity depends on the strength of the ex- 
change coupling J z (and hence the Luttinger interaction 
parameter) . 



III. QUENCH OF A SPATIALLY VARYING 
MAGNETIC FIELD IN A XXZ CHAIN: 
BOSONIZATION APPROACH 



In this section we will study a general XXZ chain which 
is subjected to a time-dependent spatially varying mag- 
netic field. The most convenient way to study this prob- 
lem is via the bosonization approach. We first review 
some of the notation. 
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FIG. 5: (Color online) Long time behavior of C xx (j = 0, n, t) 
for different values of the domain wall height mo defined as 
Sj (t = 0) \ — —sgn(j)mQ. The solid line is the equilibrium 



correlation function C%%. 



A. Equilibrium correlations 

The hamiltonian for the XXZ-sjim chain in a magnetic 
field is 



H 



(37) 

The above can be mapped onto a Luttinger liquid with 
a back-scattering potential^, 



H=— dx 
2tt 



K(V0) 2 + ^(V^) 2 



dxh(x)V<j) 



2g 



dxcos(4(/>(x)), (38) 



where g — J z a, and the fields and <f> are defined in 
terms of the boson creation/annihilation operators as 



[x)=-{N R + N L ) 



TlX 



1/2 



( X ) = (N R -N L )^ 



+ fE(^) 1/2 T^^ N/2 ^(^-M-(40) 

where N a = RiL = J2 P ^o '■ V4-(p)Vv(p) : , and : ■•• : de- 
notes normal ordering of operators. In the continuum 
limit, we interpret x = ja for some integer j and lattice 



spacing a. The coupling-parameter relationships are 



uK = vp = Jasin^kpa) 
u 



K 



1 + — — (1 - cos(2fc F a)) 
ttvf 



(41) 
(42) 



For an xy anti-ferromagnet ( J > 0), the mapping from 
spin language to boson language is given by^ 



S z {x) = — V<f>(x) + LJZcob[2^(x)], (43) 
7r 7ra 



i 



V2 



7TO 



: exp[-»0(z)] [(-If + cos[20(a?)]](44) 



In equilibrium and zero magnetic fields, the magnetiza- 
tion and basic physical correlations in the gapless anti- 
ferromagnetic XX phase are 



S'(x)) = ~(V<Kx))=0, 



S z (x + n)S z {x) 



r'j ( -) +c 2 (-i) n (- 

1 n I \n 



2 A 



S + (x + n)S-(x)) = 



C 3 - 



^ \ 2K ~< 27T 



1 \ 2K 

+ C7 4 (-1) W ( - 



(45) 



(46) 



(47) 



where the Cj are non-universal constants. 

The antiferromagnetic Hamiltonian may be trans- 
formed to a ferromagnetic Hamiltonian by the transfor- 
mation Sj — > (—iySj, resulting in ferromagnetic oper- 
ators 



S z (x) { . m . = — V#a 



-If 



cos[2<j>(x)], 



(48) 



S + {x) Lm . = 



\2ttol 



exp[-z0(x)] [1 + (-l) x cos[2<?!.(x)] 



(49) 



As a consequence, the magnetization and spin correla- 
tions in the gapless ferromagnetic XX phase are, 



1 



S'{x)) = — <V#0> = 0, 

/ f.m. T 

S z (x + n)S z (x) y 



f.m. 



2 A" 



S+(x + n)S-(x) 



(50) 



(51) 



(52) 
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B. Magnetic field quench in the XX model (K = 1) 

1. Diagonalization 

In this section we will study the time-dependent Hamil- 
tonian 



H = Hid(-t) + H f 0(t) 



(53) 



We choose a point in the XX phase corresponding to the 
XX model K = 1 and u = v F . Thus, 

Hi = — [ dxf (V6>) 2 + (V</>) 2 + — hCx)V(f> |(54) 
2tt J \ v F J 

= ^2v F \p\a\a pi (55) 

so that the initial state at t < is the spatially inhomo- 
geneous ground state of H^. The magnetic field h(x) is 
suddenly switched off at t — so that the system evolves 
according to 



Hf = v ij ^(( V0 ) 2 + w) 2 



= ^2v F \p\blb p 

p^O 



(56) 
(57) 



Since Eq. [54] basically represents harmonic oscillators 
subjected to an electric field, it may be easily diagonal- 
ized by the shift, 

(58) 



bp — a,p + Ap. 



where, 



Ap , — ■ — j — hp , 



VF\f2-K\p\L 



(59) 



and h p is the Fourier transform of h(x). Taking h(—x) = 
— h{x) (corresponding to a magnetic field that is anti- 
symmetric in space allowing for a domain wall solution), 
this gives 

hp = —i J sin(px)h(x)dx. (60) 

For some time t after the quench, it is straightforward 
to show that 

<f>(x, t) = e lH ^<t){x, 0)e-^* = <j> a (x, t) + 6<j>(x, t), (61) 
9(x, t) = e lH ^9{x, 0)e- iH f t = 9 a {x, t) + 89(x, t), (62) 

where 

X l e -»\p\/2-i P x e top\p\t + a _ pe -™ F \p\t^ ^ (g 3 ) 



mt +\ j_ i7r f l \p\ 

p^o v 

x _}_ e -a\p\/2-ipx ^ a t e ™Hpl ; 



1/2 



\P\ 



-lVF\p\t 



) ■ (64) 



Since we will be computing expectation values with 
respect to an initial state which is the ground state of 
the o p operators, <f> a ,9 a will just return the equilibrium 
results. Thus the effect of the shift and hence the initial 
magnetic field is contained entirely in 



S(f>(x 



89(x,t) 



i f°° dp 

,t) = / — cos [pvFt] cos [px]h p , (65) 

KVF Jo P 



dp 



sm[pvFt] shi[px]h p . (66) 



ttv f J P 

For any antisymmetric h(x) the above reduces to, 
1 f°° 

54>(x, t) — / dx'h(x') [sgn[a/ + z_l 

8«F J-oc 

+sgn[x' + z + ] + sgn[.x' — z + ] + sgn[x' — Z-]] , (67) 
1 



89{x, t) = — dx'h{x') [sgn[x' + z+] 
+sgn[x' — z + ] — sgn[x' + z_] — sgn[a;' — Z-]] , (68) 

where sgn[ir] = x/|x|, and z± = x ± vpt. 

The above implies that the magnetization before the 
quench is 

S') = -k^P-) = = — h(x) (69) 

/ n ox tt ox ttvf 

Thus within the bosonization approach, the magnetiza- 
tion follows the magnetic-field profile. 



2. Magnetization and Correlations after the quench 

The magnetization at some time t after the magnetic- 
field quench is 



S' 



1 



{d x <t> a {x,t)) - -d x 8<t>(x,t) 



Ittvf 



(h(x + vpt) + h(x — vpt)) 



(70) 



As anticipated, the domain wall spreads out with speed 
vf in both directions. This result is consistent with the 
ballistic time-evolution found in section [TT] from the exact 
solution of the fermionic problem. 

Next we turn to the evaluation of the C xx correlation 
function which from the exact solution of the fermionic 
model was found to be (^Sj +n Sj\ -> cos (^p) C xx (n) 

where C xx (n) ~ was the equilibrium correlation func- 
tion, and the wavelength of the oscillation was found to 
be A = too being the height of the domain wall (or 
the magnitude of the maximal magnetization at the ends 
of the domain wall) . 

In the bosonization approach, 



S + (x + n,t)S~(x,t) 



exp[-i(89(x + n,t) - 89(x,t))]-=. (71) 
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Let us suppose 



giving us precisely the observed behavior of Eq. 



h(x) — ho tanh(a;/£) 



implying 



hp = -ih 



sinh (^J 



(72) 



(73) 



Then we find, 



5d(x + n,t)-66(x,t) 



2v F 



loe 



cosh[(x — v F t)/£\ cosh[(.T + n + v F t)/!;] 
cosh[(x + n — v F t)/£] cosh[(a; + v F t)/£] 



v F t,\v F t±x\,\v F t±(x+n)\^£ 



n + v F t\ 



^0 

2v F 

Vpt\ 



V F t\ 



"1 



(74) 



For vpt >> \x\, \x + n\, i.e, long after the moving domain 
wall front has crossed the two observation points at x 
and x + n, we find 



<J0(a: + 7i,i)-<S0(a;,*) 



Vp 



so that 



(75) 



(76) 



Thus we find that for an initial domain wall that has 
the profile 



(S z {x = ±oo)) 



± 



ho 



TTVp 



±m 



(77) 



the transverse spin correlation function at long times do 
not equilibrate but acquire a spatial oscillation at the 
wavelength 



A 



2nvp 



2 

m 



(78) 



This result is identical to that obtained from the exact 
solution of the fermionic problem (Eqns. [34l |36|) . 

Note that the wavelength of the spatial oscillation is set 
by the height of the domain wall mo and is independent 
of the width £ (or a = in the fermionic problem). 
For the case of a linearly varying magnetic field studied 
in the previous section, the domain wall magnetization 
approached its maximum possible value at its two 
ends. This in the boson language corresponds to 



ho 

Vp 



TV 

"2a' 



(79) 



so that writing n — na in Eq. 1761 where n is an integer, 

' ' (80) 



(S+(x + na)S-(x)) -> e l7rfi/2 -= 



For the evaluation of the (S Z S Z ) correlation function, 
we need to evaluate exp [iS<p(x + n, t) — i6(f>(x, t)]. We 
find, 



4>{x + n,t) — 64>(x, t) 
-hoi 



2v F 



loe 



cosh 


e 


cosh 


x+n— v F t 


cosh 


X-\-Vpt 

s 


cosh 


X — Vpt 





(81) 



v F t^>£,\x\,\x+n\ 



> 



The above result together with the fact that 



d x S(f>(x + n, t)d x S4>(x, t) 



A-kv f 



(h(x + vpt) + h(x — vpt)) 



x (h(x + n + vpt) + h(x + n- v F t)) . (82) 



implies that for v F t ^> £, \x\, \x + n\, the C zz correla- 
tion function equilibrates which is in agreement with the 
results in section HTl 

To summarize the results of this section: a). Within 
the bosonization approach, an applied magnetic-field pro- 
duces a spin pattern that follows the local magnetic field, 
b). After a quench, the domain wall spreads out bal- 
listically. The magnetization equilibrates (in this case 
vanishes to zero everywhere), c). The C zz correlation 
function equilibrates, whereas the C xx correlation func- 
tion acquires spatial oscillations at wavelength A = — 
where mo is the height of the initial domain wall defined 
as iSz ( ± oo ) = ±mo- These results are in agreement with 
the exact results of section [TT] obtained from solving the 
fermion lattice problem. Thus bosonization is good at 
capturing the behavior at long distances and times, even 
for the nonequilibrium problem. As expected it misses 
some of the details of the short distance physics both 
in the static properties (as shown in Fig. [6]), as well as 
during the time-evolution such as the existence of a com- 
plex internal structure in the propagating domain wall 
fronts. 

In the next section we provide an explanation for the 
inhomogeneities observed in the transverse spin correla- 
tion function by solving the classical equations of motion 
for the spins. 
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FIG. 6: (Color online) Plot of the magnetization obtained 
from the bosonization approach with h(x) — ho tanha;/i;, and 
from the fermionic lattice model with a — 20. 



FIG. 7: (Color online) Thick line: Time evolution of the mag- 
netization from the solution of the classical Landau-Lifshitz 
equation assuming that S z (t < 0) is given by Eq. Thin 
lines: Time evolution from the exact solution of the quantum 
problem. 



IV. MAGNETIC FIELD QUENCH IN A 
CLASSICAL XX SPIN CHAIN: FORMATION OF 
A TRANSVERSE SPIN- WAVE PATTERN 

The Hamiltonian in Eq. [1] implies the following equa- 
tions of motion for the spins 



dS 



dsi 

dt 

ds] 



-JS 



3 (s y 



J'+l 



J'-l/ 



hjSf 



(83) 
(84) 

(85) 



We will treat the spins as classical variables, and solve 
the above equations for a sudden quench of a spatially 



inhomogeneous magnetic field, hj(t) 



-t). First 



notice that for t < 0, after performing the course-graining 
Sj = Sj +1 , the following spin profile satisfy the equations 
of motion, 



b > ' 2J 



S* = 



S 2 - (S 



z )2 

j I 



s» = o 



(86) 

(87) 
(88) 



With the above as an initial condition, the dynamics 
for t > can be easily studied numerically and we find 
a ballistic spreading out of the domain wall profile (see 
Fig. [7]), and the appearance of a spin- wave pattern. The 
latter result, i.e., the appearance of an inhomogeneous 
pattern at long times can be verified analytically rather 



easily as follows. Suppose, hj(t) 



-ho tanh | 



9(-t). 



Then for t > 0, if we assume a ballistic broadening of 
the domain wall, 



SUt) = 



hp 
2J 



, . ja — vt. , ,ja + vt. 

tanh(- — ) + tanh( : ^ — - — ) 



(89) 



Writing Eq. EH [85] in terms of S ± = S x ± iS y , and 
performing a course-graining, one obtains the following 
equation for S + (x = ja) 



dS+(x,t) 
dt 



2iJS z (x,t)S + (x,t) 



(90) 



Using Eq. [89l Eq. [90] may be easily solved to obtain 



S + (x,t) = 5+(a;,0)exp 



ih Q £ (, cosh ( 



vt-\-x 



2v 



In- 



cosh(- 



(91) 



Thus for vt 3> | a; | we find the spin- wave pattern 

S + (x,t) -> S + (x,0y hox/v (92) 

in agreement with the results of the previous sections. 

Thus the appearance of the spatial inhomogeneity may 
be understood as follows. Right after the magnetic field 
is switched off, the spins in the XY plane begin to precess 
due to a non-zero magnetization in the z-direction. How- 
ever this precession only lasts for as long as it takes the 
domain wall profile to flatten out to zero. Since this time 
is different for spins located at different spatial positions, 
the spins locally dephase with respect to each other and 
arrange themselves in a spin-wave pattern. 
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V. QUENCH OF A SPATIALLY VARYING 
MAGNETIC FIELD AND INTERACTIONS IN A 
XXZ CHAIN: BOSONIZATION APPROACH 

A. Quench within the gapless XX phase 

We will now turn our attention to the case where as 
the magnetic field is switched off at t = 0, the interac- 
tions are also turned on so that the Luttinger parameter 
is K — 8(—i) + K6(t). We will first consider the case 
where the ground state both before and after the quench 
corresponds to a gapless XX phase, so that the cos(/3</>) 
term can be neglected. 



The §(</>, 0) arise due to the spatially varying magnetic 
field and are given by, 

u 

5c/>a(x,t) = — 6</>(x,t,P = 0,v F -t u), (103) 

Vp 

ue 20 

60p(x,t) = 60(x,t,/3 = 0,v F -> u), (104) 

Vf 

where 8cj){x,t,j3 = 0, Dp), S0(x, t, /3 = 0,vp) are given in 
Eq. [53] and [55J respectively. 



2. Magnetization after the quench 



oliubov rotation 



We wish to begin with a non-interacting system in an 
inhomogeneous magnetic field, 

Hi = £ / dx ( (v< ^ )2 + (w)2 + ^- /i ( a; ) v ^)( 93 ) 

= ^2v F \p\a\a p , (94) 

and quench to zero field, while turning on interactions so 
that 

Hf = lr I dX (^ (V0)2 + K(V9)2 ) (95) 



(96) 



The rip operators have a simple time evolution and may 
be related to the b p operators by a Bogoliubov rota- 
tion^, 



i] p = cosh (3b p + sinh/36L p , 
if v = cosh + sum /36_p. 

where e~ 2 ^ — K. We define, 



(97) 
(98) 



f(p,t) = cos(u\p\t) -ism(u\p\t) cosh(2/3), (99) 
g(p,t) = isin(M|p|t)sinh(2^). (100) 

and a± = f(p,t) ± g{p 1 t). Performing the shift b p — 
a p + A p as before, we write the 4> and 9 fields as (cj), 9) = 
(tt%,0%) + 5{<t>f,,Op), where 



,a, .s V"^ / ^P 



L V 2tt I p 



_ -i7T^{L\p\Y /2 1 -a\p\/ 



12-ipx I* J 



+a + a_ p ) 



(101) 



p^Q V 7 U ' 



(102) 



It is straightforward to show that the magnetization at 
a time t after the quench is 

S z {x,t) = = _L_ + „t) + - u t)] 



7r dx 2ttvf 



(105) 



Thus even with the interaction quench, the domain wall 
spreads out ballistically, but with the renormalized ve- 
locity u = vf/K. 



3. Basic non- equilibrium correlations 

To compute the correlation functions, some basic ex- 
pressions we need are 



G^ = <[^(x J t)-^(x / ,t)] a > 
Gee = {[0}{x,t) - 0}{x\t)] 2 ) 



(106) 
(107) 



Note that Gee is related to G^ by K — > -k the com- 
putation of which is identical to that already presented 
i n 17 ' 18 . For completeness we present the results below, 



G 00 = / — [cos 2 (wtp) + e~ 4/3 sin 2 (utp)] 
Jo V 



x(l-cos[p(x-x')})e~ ap 
which eventually gives 



(108) 



8 



loe 



[a 2 + (2ut) 2 ] 2 [a 2 + x 2 ^ 2 



Gee — G 



loe 



a 4 [a 2 + (2ut + x) 2 ][a 2 + (x - 2ut) 2 ) 
K- 2 -Y 



(o) 
ee 



[a 2 + (2ut) 2 } 2 [a 2 + x 2 } 2 



_a 4 [a 2 + (2ut + x) 2 ] [a 2 + (x - 2utf 



(109) 



(110) 



where G 



(0) _ K 



log 



,G 



(o) _ j_ 

IK 



loe 
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4- Spin correlations after the quench 

Taking the 2ut » x limit, we can easily obtain the 
correlator asymptotics. The memory of the initial do- 
main wall profile appears as oscillatory pre-factors of 
the form exp[2i(S(j)(x + n) — S<f>(x))] and terms such as 
V8<t>(x + n)V5(f>(x). For the specific case of h(x) = 
ho tanh(x/£) we find at long times, 



(111) 



S z (x + n)S z (x)) = f(x + n, t)f(x, t) 

r,2 



,iC 2 + l 



S+(x + n)S-(x) 
Ci(-l) n e^ 



A+K 2 + 



l + K- 2 



where 



f(x,t) 



2v F 



(h(x + ut) + h(x — ut)) 



(112) 



(113) 



The case of the pure interaction quench (ho — 0) was 
studied i n 17 ' 18 where the above anomalous power-laws 
corresponding to an algebraic decay which is faster than 
in equilibrium (since K 2 + 1 > 2K for all real K) were 
obtained. From Eq. 11121 we find that the effect of an 
initial domain wall profile is to superimpose spatial os- 
cillations onto this decay, where the wavelength of the 
oscillations depends on the interaction parameter K and 
is found to be 



A 



2-kvfK 2K 



ho 



mo 



(114) 



As expected, increasing the exchange interaction J z in- 
creases the effective magnetic field seen by the precessing 
spins, thus decreasing the wave-length of the spin wave 
pattern. 



5. Generalized Gibbs Ensemble 

Rigol et al have made the interesting proposal that the 
long-time properties of out of equilibrium integrable sys- 
tems may be captured by a generalized Gibbs ensemble 
that enforces the constraints imposed by the integrals 
of motion—. This approach has been applied with suc- 
cess by Iucci and Cazalill a 17 ' 18 who studied homogeneous 
interaction quenches in models similar to those studied 
in this paper. It has also been successfully applied to 
models of free bosons 3 ^., and for some general family of 
integrable models^ 7 -. An interesting question has been 
to identify situations in which the Gibbs ensemble ar- 
gument might break dow n 38 i 39 . We find that for an in- 
homogeneous quantum quench starting from an initial 
domain wall state, the Gibbs ensemble argument cannot 



capture the spatial inhomogeneities at long times after 
the quench. In this section we review the Gibbs ensem- 
ble argument and show why this does not work for our 
case. 

Consider a quantum quench of the form H = 
Hi6(—t) + Hf9(t), so that for t < the system is in the 
ground state of Hi, whereas at t > 0, the wave- function 
begins to evolve according to Hf. Let us further assume 
that both Hi and Hf can be diagonalized 



Hi = V elala. 



H 



f 



p 



(115) 
(116) 



b p = cosh Opdp 



and that the two Hilbert spaces may be related to each 
other via a canonical transformation. For the case of 
bosons a general canonical transformation is of the form, 

(117) 

where X p is a linear shift and p is a rotation angle . 

The essential idea behind predicting the long-time be- 
havior is that an integrable model such as Hf has a con- 
served quantity for every degree of freedom. For the ex- 
ample given above, the conserved quantity is trivially 
given by I p — b p b p . Thus during the time-evolution 
(bpbp) should be conserved. Since the initial state is the 
ground state of Hi for which (aj,a p ) =0 for p ^ 0, it 
follows that the Gibbs' ensemble should be characterized 
by the distribution function 

(118) 



(blbp) = |A p | 2 +sinh 2 e f 



For the particular models studied here, we find the 
initial spatial inhomogeneity is captured by the linear 
shift X p (given by Eq. [59}, while the angle p captures 
the homogeneous changes in the interaction parameter. 
Since in our example |A P | 2 vanishes in the thermody- 
namic limit, any effect of the initial spatial inhomogeneity 
is lost. Thus while the Gibbs ensemble captures the ef- 
fect of the rotation, and therefore correctly predicts the 
new power law exponents o£i£, it misses the spin-wave 
pattern at long times. In particular the Gibbs ensemble 
prediction for the transverse spin correlation function is 
given by 



S + (x + n)S-(x) 

cu-i) n 



c 3 



(119) 



which differs from the correct answer (Eq. I112j) by the 
absence of the spatially oscillating factors e v f K . 



B. Magnetic field and Interaction quench into the 
gapped Ising phase: Semiclassical treatment 

In this section we will study the effect of quenching 
the external magnetic field and the interactions into the 
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massive Ising phase where the cos <f> term is strongly rele- 
vant. In order to study the time evolution, we will employ 
a semi-classical approximation which corresponds to rc- 
placing cos0 ~ 1 — Thus for this case the initial 
wave-function is the ground state of 



Hi 



v F 



dx 



2tt 



(v^) 2 + (w) 2 



V F 



-h(x)V(j) 



(120) 
(121) 



p 



while for t > the wave-function evolves according to 

™2 



. dx 

2tt 



(v</>) 2 + (yey 



(122) 
(123) 



where (j> = -J= , 9 = \f~K9 and co p — yjp 2 u 2 + m 2 . 

The calculations mirror that of the previous section. 
Writing (f> = (f> a + 8(f), 9 = 9 a + 89, we find for a time t 
after the quench 



p^O 



1/2 



x _ e -<x\p\/2-ipx 
P 



a p ^cos(ojpi) + iK sm(uj p t)j 



Mp\ 



+a_ p I cos(wpt) — iK sin(a;pt) 



(124) 



1/2 



-a|p|/2— ipa; 



a| ^cos(wpt) + s j^jpj sin(wpt) 



-a_„ ( cos(wpt) — i ^— - sin(a;pt) 

' 1 ' Au|p| 



(125) 



and 



% f°° dp 

84>{x,t) = / — cos [uipt] cos [px]h v , (126) 

irv F Jo P ' 

S9(x,t) = — - [ — sin[w„t] sin[px] f ] h p 
nKv F J p \u\p\J 

(127) 



1. Magnetization after the quench 

The magnetization at a time t after the quench is given 

by (S z (x,t)) = -id^OM) + ^(cos^OM))}- Due 
to the presence of infra-red divergences, we find that for 
all times after the quench, 



so that Ising order never develops. This is consistent 
with the results of-2-2. and may be an artifact of the semi- 
classical approximation. In section I VII we will compare 
this result with a mean-field treatment for the order- 
parameter dynamics. 

Thus, within the validity of the semiclassical approxi- 
mation we find the following time evolution for the initial 
domain wall profile, 



(S*(x,t)) 
ho 



1 



7T 



d x (8<f>(x,t)) 



dz cos I — • I ^ 2 

KVF Jo \ £ 



z z + 



u 



sin 



O) 



sinh(f) 



(129) 



The above is plotted in Fig.|S]for several different masses. 
The domain wall profile shows a more complex time evo- 
lution than for the quench in the gapless phase where a 
simple ballistic broadening of the domain wall occurs. In 
the Ising phase instead we find oscillations of the local 
magnetization at the frequency of the mass. In addition 
the domain wall spreads out less and less as the mass 
increases, a result which is in agreement with a t-DMRG 
study of the XXZ chain 12,13 . Note however that we find 
oscillations of the magnetization outside the light-cone at 
the time-scale of 1/m. This is due to the nature of the 
Klein-Gordon equation d\ — <9 2 + to 2 which has solutions 
of the form e ikx-Wk»+mH wt>l ) e _ im4 jkx-if^t _ The t _ 
DMRG results on the other hand found oscillations only 
within the light cone. 

We also find that the magnetization eventually equi- 
librates at long times. Fig. [S] shows the time evolution 
for the magnetization for a given point in position space. 



At long times Eq. 11291 may be evaluated in a stationary 
phase approximation which gives, 



^(■'■")-^V^ 00B ( m * + i) (130) 



itvf 



showing that the magnetization decays locally as a 
power-law. 



2. Correlation functions after the quench 

Here we give results for some typical spin correlation 
functions. Again, due to the presence of infra-red diver- 
gences, the transverse-spin correlation function is found 
to vanish at all times after the quench, 



I e i9(x,t) e -i6{0,t)x ^ q 



(131) 



On the other hand the z — z correlation function is found 
to be 



(128) ^ e 2i0(x,t) e -2#(O,t)^ = e M(S«x,t)-6K*,0)) Al ( a ) Q^132) 



14 



m^/u = 




20 30 
time ut/^ 



FIG. 8: Contour plot for \S Z \ after a "semi-classical" quench 
for several different masses. 
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FIG. 9: (Color online) The magnetization decays as t 1 ' 2 
with oscillations of period ~ m _ . Here, m£/u — 0.1. 



For a purely interaction quench (50 = 0) we find a power 
law decay in the Ising phase with an exponent which is 
half of that in the initial gapless phase, a result which was 
also found by-i& for homogeneous quenches. The initial 
domain profile superimposes the additional structure 



a 2i(S<t>(_x,t)-S<f>(y,0)) 2Hl 



(1133) 



Thus at long times, the spatial oscillations due to the 
initial domain- wall eventually decay to zero. 



VI. EFFECT OF QUENCHING A MAGNETIC 
FIELD AND INTERACTIONS: MEAN-FIELD 
TREATMENT 

In this section we will use a mean-field argument to 
explore whether after quenching both a magnetic field as 
well as interactions, local anti-ferromagnetic (AF) Ising 
order can develop. Thus our initial wave-function corre- 
sponds to the ground state of the following Hamiltonian 



QX QX , QV QV , "3 q2 



(134) 



with J > 0, and we will assume hj = Tja. While, as 
a result of the quench the wave-function evolves in time 
due to 

h, = jJ2 (s*s* +1 + S y S y +1 ) + j'Y, s-s> +1 (135) 



Performing a Jordan- Wigner transformation followed by 
the change of variables Cj — > (— l) l Ci, Hi becomes iden- 
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tical to the Hamiltonian (Eq. [2J already studied in sec- 
tion whereas the Hamiltonian after the quench is, 



J 



^ (c\c i+ i + h.c. 



(136) 



Thus the initial state which is the ground state of Hi 
(given by Eq. [7]) evolves in time according to Eq. 11361 
In what follows we will study this time evolution using a 
mean-field approximation. 

Defining the AF order-parameter at a time t, Aq (t) as 
(4(t)ci(t)) = ± + (-l) l A (t), Eq.HSf)] within a mean-field 
approximation becomes 



H 



mf 



£k ( C k Ck - 4+7T C fc+7r 



| fc| <tt/2 



-2J*A (i) Y, (4 



"/C+7T 



Cfc 



(137) 



| fc| <tt/2 



where e k = — Jcoska, and Ao(£) is determined from the 
self-consistency condition 

Ao(t) = ^me^f'Yi-iyche-^T^) (138) 

i 

where N is the total number of sites. 

In terms of the amplitudes c k ± of being in the adia- 
batic eigenstates of H J 1 , the self-consistency condition 
becomes 

A °W = ^ E [^29 k (t)(\c k+ (t)\ 2 -\c k .{t)\ 2 ) 

|fc|<7r/2 

-cos29 k {t) (<? k+ {t)c k -{t) + c* k _(t)c k+ (t)) x 
cos( / dx{E k+ (x) -E k _(x))) 



(139) 

(140) 
(141) 



-icos29 k (t) (c* k+ (t)c k -(t) - ct_(t)c k+ (t)) x 



sin( / dx(E k+ (x) - E k _(x))) 
Jo 



where 



E k± (t)=±yJe 2 k + U2A 2 (t) 
2J z A (i) 



tan26> fc (i) 



e-k 



and the coefficients c k ± satisfy the equations of motion, 



dc k ± (t) 
dt 



= c kT (t) 



±e k J; 



dA„ 
dt 



4J?A2 



,±if*dx(E h+ (x)-E h _(x)) 



(142) 



In an adiabatic approximation A§ 1, we can ap 
proximate Cfc± by its value at t = 0. In this case Eq. [T 



becomes 

A ° " E l sin2 °k(t) (sm29 k {t)I£ + cos28 k (t)r k ) 

| fc| <tt/2 

-cos26> fe (t) (sin20 fc (t)/| - cos 26> fc (£)/£) x 



cos( / dx(E k+ (x) - E k -{x))) 



-icos26» fc (i)^sin( / dx{E k+ {x) - E k _(x))) 
Jo 



(143) 



where I k ' y ' denote the following expectation values in 
the initial state before the quench, 



i k = (*l4 c fc+- + 4+^110 



(*|CfcCfc 



-"fc + 7T 



Cfc+Trl*) 



(144) 
(145) 
(146) 



Since in the initial domain wall state, I^' y ' z ~ 0, this 
implies that Ao(i) = 0, and hence anti-ferromagnetic or- 
der does not develop within a mean-field and adiabatic 
approximation. This result can also be understood from 
the fact that the initial state is a highly excited state of 
the anti-ferromagnetic model as it corresponds to intro- 
ducing a domain wall in the staggered magnetization at 
almost every sitei£. 

The above conclusion also holds if we relax the adi- 
abatic approximation. To see this note that at every 
fc, the mean-field Hamiltonian Eq. 11371 corresponds to 
a pseudo-spin 1/2 r k where^ t\. 



c k c k C fc+7r Cfe+ w , 



4 + 71-' 



rl.^Ck and T y 



(4+^ 



This pseudo-spin is subjected to a magnetic field that 
has to be determined self-consistently at each time t 
where for t = (r*' z (0)) = I%' z , <r« (0)) = iJ* Since 
for the initial domain wall state (r(0)) = 0, it implies 
d }jn^ (t = 0) = to all orders in n, so that anti- 
ferromagnetic order (oc J2 k ( T k)) wm " n °t develop. 



VII. CONCLUSIONS 

In summary we have studied quenched dynamics in a 
XXZ chain starting from an initial inhomogeneous state 
where the spins are arranged in a domain wall profile. 
The dynamics of the domain wall is found to be quali- 
tatively different within the gapless XX phase and the 
gapped Ising phase. In the former the domain wall 
broadens ballistically, while in the latter the domain wall 
spreads out less and less the deeper one is in the Ising 
phase. The magnetization is locally found to oscillate, 
and eventually decays as a power-law. Although the re- 
sults in the Ising phase were obtained within a semiclas- 
sical treatment, they are in qualitative agreement with 
t-DMRG simulations of the XXZ chain^^. 

We have also presented results for the time-evolution of 
two-point correlation functions. For quenches within the 
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gapless XX phase we find that the longitudinal spin cor- 
relation function (C zz ) reaches a steady state which is in- 
dependent of the initial domain wall. On the other hand, 
the transverse spin correlation function (C xx ) reaches a 
nonequilibrium steady state which retains memory of the 
initial domain wall. This memory appears as an inho- 
mogeneous spin wave pattern with a wavelength which 
is inversely related to the height of the initial domain 
wall. This result was obtained by using three differ- 
ent methods: exact solution of the fermionic problem 
(Eqns. El [Ml), bosonization (Eqns ITT21 [TTIj) and the 
solution of the classical equations of motions for spins 
(Eq. [92]). The spatial oscillations are found to arise due 
to the dephasing of transverse spin components as the 
domain wall broadens. We also find that a Gibbs ensem- 
ble argument is not adequate in capturing the spin wave 
pattern. 

For quenches into the gapped Ising phase, all inho- 
mogeneities both in the local magnetization and in the 
two-point correlation functions are found to eventually 
decay away. This may be an artifact of the semiclassi- 
cal approximation which neglects soliton creation^ We 
would also like to point out a recent preprint — that stud- 
ies time-evolution of an initial domain wall in the gapped 
Ising phase using Algebraic Bethe Ansatz. The authors 
study the Loschmidt echo and find a lack of thermaliza- 
tion. 

Finally we perform a mean-field treatment for the 
time-evolution of the anti-ferromagnetic Ising gap and 
find that for an initial state corresponding to a do- 
main wall profile, anti-ferromagnetic order never devel- 
ops. This should be expected on the grounds that the 
initial domain wall profile corresponds to a highly excited 
state of the anti-ferromagnetic Ising model. 
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tors S • = ±iSj as follows, 



S'7 



exp 



exp 



-«r E c «' 

n=-N/2+l 



E 4< 



n=- N/2+1 



(A2) 



(A3) 



The Hamiltonian is diagonal in momentum space Cj = 
•JR Efc c fce lfcj where, 



H xx = ^ £kc{ 

k 

€k = — J COS k. 



Ck 



(A4) 
(A5) 



Thus the operators Cfc, have the trivial time-evolution 
Cfc(t) = Ck(0)e~ lekt . In terms of the r] m quasi-particlcs 
that diagonalize the Wanner-Stark problem (Eq. [5]), the 
transformation given by Eqns. (0 E} yields 



k,m 



e ik(j- m )-ia sin k e -ie k t Vm ( A g) 



e -ik(j-m)+iasuik e ie k t T j\ (A7) 



Appendix B: Time evolution in the XX model 
starting from initial domain walls of different heights 

Antal et al. have considered more general domain 
walls 7 in which the homogeneous magnetization on either 
side of the wall is taken to be ±mo, where < too < 1/2. 
They have studied the time evolution of the domain wall 
under the influence of the XX hamiltonian. In this sec- 
tion we will extend their results by studying how the 
transverse spin correlation function evolves in time. 

A domain wall of height too may be constructed as 
follows 7 -, 



l*mo>= 114 11 4 10} 



(Bl) 



k= — k— k— — k^ 



Appendix A: Diagonalization of the XX-model 



The Hamiltonian for the XX model is 



H , 



3 



(Al) 



The Jordan- Wigner transformation maps the spin Hamil- 
tonian to that for spin-less fermions Cj and c| which are 
defined in terms of the spin raising and lowering opera- 



wherc 



■3' 



j>0 



3 ' 



3<0 



(B2) 
(B3) 



and k± = it ± mo) . Due to the XX Hamiltonian H = 
— J^2 k cos(k)cj,Ck, the operators Cj evolve as 



i (t)=J2i m ^J m _ j (Jt)c m (Q), 



(B4) 
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In order to compute any correlation function, the con- 
traction we need is 



J.(t) Cj+n (t)) = J2 iS ~ m ~ nj ™-j{Jt)Js-3-n{Jt) 

x (4(0)0,(0)). (B5) 



where, 



if s ■ m < 

, r t (n , iJ(l))H ^Ki±p^)] if s , TO<0 
sin[(f ffmo)(8 - m)1 

7r(s — mj 



Thus, 



(B6) 



c](t)c j+n (t)) = 

711 , S < 

sin [(| + 7rm ) (m — s)] 
ir(m — s) 



m,s>0 

sin [(J — 7rmo) (m — s)l 



(B7) 



7r(m — s) 

We write both terms as sums over positive to, s, (i.e., 



letting m, s — > — to, — s in the second term), and rewrite 
m>o ES-oo with I — s ~ m, and use the 

identity 



>/j+m(^)>/j+n+m(-T) — — (J j+l(x) J j +n (x) 

m>0 

-Jj(x)J j+n+l (x)). (B8) 



Next we employ the asymptotic expansion for large ar- 
guments, J n (x) ~ \Jljjj: cos [ x — f n — f ] to nn( l the fol- 
lowing nonequilibrium steady state expressions, 

lim (cUt)c j+n (t)) ~ 

-n Sin [(f ~ 7rm o) '] Sin [f ~ n )] 



E 



ir 2 l{l-n) 



I — — OO 

■i+n sin [(f + tttoq) I] sin [f (j + n)] 



tt 2 Z(Z + n) 



(B9) 



Examining the limit of too = i, we recover exactly the 
values in Eq. [28l [29] We use Eq. IB9I in the evaluation 
of the determinant in Eq. [23] required for computing the 
transverse spin correlation function. We recover a spatial 
period of oscillation in C xx (j, j +n, t — > oo) that increases 
with smaller too (see Fig. [5]) . 
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